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Preface 


This  study  was  the  first  attempt  at  providing  a usable  set  of 
computer  routines  (which  employed  current  concepts  for  the  selection 
of  the  design  space,  second  order  surface  fit  approximations,  and  a 
conjugate  gradient  minimisation  algorithm)  to  be  used  by  the  Air  Force 
Flight  Dynamics  Laboratory  Design  Branch  in  the  initial  design  phase  of 
an  aircraft.  Several  proposals  and  programs  were  available  but  the 
lack  of  information  as  to  the  use  and  application  of  these  concepts 
led  to  this  study.  Computer  routines  that  had  been  working  were  found 
not  to  be  compatible  with  the  Wright-Patterson  Air  Force  Base  computer 
system  and  had  to  be  modified  and  rewritten.  This  study  attempted  to 
eliminate  these  problems  and  to  demonstrate  the  usefulness  of  the 
surface  fit  approximation  scheme. 

I would  like  to  thank  Mr.  Jim  Parker  and  Captain  Russ  Morrison  of 
the  Air  Force  Flight  Dynamics  Laboratory  Design  Branch  for  their  assist- 
ance in  providing  me  with  input  Information  for  the  mission  simulation 
used  in  this  study. 

1 would  also  like  to  thank  Dr.  K.  S.  Nagaraja  of  the  Air  Force 
Flight  Dynamics  Laboratory  Aerodynamics  Branch  and  Major  Steve  Koob  of 
the  Aeromechanlcal  Engineering  Department  of  the  Air  Force  Institute 
of  Technology  who  provided  the  thesis  topic  and  valuable  aid  during 
the  course  of  this  study. 

Finally,  I would  like  to  express  my  appreciation  to  Major  Gerald 
Anderson  of  the  Mechanics  Department  of  the  Air  Force  Institute  of 
Technology  who  agreed  to  sponsor  this  thesis  topic. 
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Abatract 

A study  was  mads  to  demonstrate  the  feasibility  of  using  surface 
fit  approximations  In  the  design  analysis  of  an  advanced  tactical 
fighter.  Design  variables  were  selected  and  a design  space  was  defined 
based  on  a simple  latin  square  method.  The  take-of f-gross-welght  (TOGW) 
and  landing  distance  (D1STL)  for  each  configuration  were  determined  by 
use  of  a computer  program  which  simulated  the  required  mission.  A 
regression  analysis  was  performed  on  this  data  to  obtain  a quadratic 
surface  fit  representing  TOGW  and  D1STL  as  functions  of  the  design 
variables.  A conjugate  gradient  technique  was  employed  to  find  the 
minimum  TOGW  subject  to  constraints  on  DISTL  within  the  design  space. 
This  value  was  compared  to  the  minimum  obtained  by  running  all  possible 
combinations  of  design  variables  through  the  mission  simulator.  It  was 
concluded  that  the  surface  fit  approximation  based  on  the  latin  square 
design  space  did  yield  accurate  results  and  required  fewer  computer 
simulations. 
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A STUDY  OF  THE  METHOD  OF 
SURFACE  FIT  APPROXIMATIONS 
FOR  AIRCRAFT  OPTIMAL  DESIGN 

I.  Introduction 

The  Air  Force  Flight  Dynamics  Laboratory  (AFFDL)  Design  Branch 
became  Interested  In  several  proposals  concerning  the  selection  of  an 
airplane  engine  by  optimization  of  surface  fit  approximations.  The 
original  studies  had  been  conducted  by  the  Aero  Propulsion  Laboratory 
(APL)  at  Wright-Patterson  Air  Force  Base.  The  results  of  these  studies 
shoved  that  considerable  savings  could  be  made  in  engineering  and  com- 
puter time  by  using  limited  quantities  of  simulation  data  (Ref  7.599). 
The  AFFDL  Design  Branch  was  Involved  in  a design  scheme  similar  to  those 
used  in  the  aforementioned  studies.  The  Design  Branch  wanted  to  inves- 
tigate  the  use  of  surface  fit  approximations  based  on  the  latin  square 
design  selector  and  Integrate  these  Ideas  into  their  design  analysis. 

In  particular  the  Design  Branch  wanted  to  apply  the  method  to  an 
aircraft  that  would  be  capable  of  accomplishing  an  interdiction  mission 
that  might  be  required  of  an  advanced  tactical  fighter.  The  vehicle 
was  to  have  a one  man  crew,  carry  5500  lbs  of  stores,  and  be  capable 
of  a supersonic  range  of  up  to  350  nautical  miles. 

The  purpose  of  this  study  was  to  demonstrate  the  feasibility  of 
using  surface  fit  approximations  in  the  design  analysis  of  an  advanced 
tactical  fighter.  The  design  process  that  was  considered  consisted  of 
six  major  steps.  The  first  step  required  the  selection  of  independent 
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design  variables  (engine,  airframe,  and/or  mission)  and  dependent  per- 
formance  or  response  functions.  The  performance  functions  would  be 
calculated  based  on  the  Independent  design  variables.  The  independent 
variables  selected  were  the  wing  loading  (airframe  variable),  aspect 
ratio  (airframe  variable),  vehicle  thrust  to  weight  ratio  (engine  vari- 
able),  and  the  supersonic  dash  range  (mission  variable).  The  dependent 
performance  functions  were  the  take-off -gross  weight  and  the  landing 
distance.  The  next  step  was  to  choose  values  for  the  Independent 
variables.  This  was  done  by  establishing  the  upper  and  lower  limits  on 
each  variable  and  then  dividing  that  range  into  an  equal  number  of 
parts  (equal  to  the  number  of  independent  variables).  The  ranges  of  the 
independent  variables  were  chosen  so  that  the  performance  functions 
could  be  adequately  represented  by  quadratic  equations  In  the  regions 
of  interest.  The  third  step  was  to  select  from  all  possible  combina- 
tions of  the  independent  parameters  that  subset  which  would  still  be 
representative  of  the  complete  design  space.  This  reduces  the  number 
of  data  points  needed  and  also  reduces  the  computational  time  for 
deriving  the  performance  function  values.  This  was  accomplished  by  a 
statistical  technique  called  the  simple  latin  square  design  selector. 

This  technique  joins  together  permutations  of  the  various  values  of  the 
independent  variables  in  such  a way  that  each  variable  is  stepped  through 
all  its  values  once  every  "1”  data  points  where  "ln  is  the  number  of 
values  for  each  variable.  The  fourth  step  was  to  compute  the  dependent 
performance  parameters  by  lnputlng  the  selected  design  variables  into 
a mission  simulator  and  outputing  the  performance  values.  The  available 
mission  simulator  was  that  used  by  the  AFFDL  Design  Branch  (the 
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Computerised  Initial  Slse  Estimate  - CISE).  This  program  is  used  to 
give  Initial  values  of  the  design  parameters  at  the  "sero  level"  of 
analysis,  i.e.,  even  before  a drawing  of  the  configuration  Is  proposed 
(Ref  4*2).  The  program  performs  a weight  oriented  siting  process  to 
produce  the  basic  physical  characteristics  of  the  vehicle.  The  next 
step  was  to  perform  a regression  analysis  on  the  performance  data  to 
yield  quadratic  surface  fits  representing  the  performance  functions  In 
terms  of  the  independent  design  variables.  Finally,  a constrained 
minimisation  of  the  performance  fmctlon  of  primary  interest  was  per* 
formed  and  located  those  values  of  the  Independent  parameters  that 
yielded  the  constrained  minimum. 

The  number  of  independent  variables  In  this  study  was  taken  to  be 
four  with  each  variable  taking  on  five  values  each.  The  number  of  vari- 
ables was  limited  so  that  all  possible  design  combinations  could  be 
calculated  and  input  into  the  mission  simulation.  The  minimum  design 
point  from  the  complete  set  could  then  be  compared  to  that  found  by 
using  the  surface  fit  approximations. 


II.  Design  Space  Selection 


General 

The  first  part  of  the  study  was  to  Identify  those  variables  that 
might  be  significant  in  the  design  of  an  aircraft  given  a general 
mission.  A distinction  is  made  between  dependent  and  independent  vari- 
ables.  The  dependent  variables  are  selected  by  considering  the  total 
design  goal.  The  TOGW  is  generally  used  as  one  dependent  variable  since 
it  can  be  calculated  from  the  various  other  parameters  by  using  a mis* 
slan  simulation.  The  Independent  variables  are  chosen  based  on  the 
ability  of  the  designer  to  change  these  parameters  and  their  Impact  on 
the  total  design  objective.  These  parameters  (both  dependent  and  inde- 
pendent) are  typically  engine  variables  (bypass  ratio,  thrust  to  weight 
ratio,  etc.),  airframe  variables  (wing  loading,  aspect  ratio,  leading 
edge  sweep  angle,  etc*),  or  mission  parameters  (range,  altitude,  loiter 
ability,  etc.). 

In  this  study,  it  was  desirable  to  have  at  least  one  variable  from 
each  of  the  three  possible  sets  of  airframe,  engine,  and  mission.  The 
mission  parameter  selected  was  the  supersonic  dash  range  (RN).  It  was 
chosen  because  it  could  be  varied  in  the  CISE  program  and  is  a 
parameter  of  interest  in  the  advanced  tactical  fighter.  The  engine 
parameter  that  was  chosen  was  the  thrust  to  weight  ratio  (TW),  i.e. , 
the  ratio  of  the  Installed  engine  thrust  to  the  vehicle  weight.  Two 
airframe  parameters  were  selected.  The  wing  loading  (WOS)  and  aspect 
ratio  (AR)  were  chosen  from  this  group.  They  represent  two  character- 
istics that  are  available  for  use  in  the  CISE  program. 


4 


r 

I 

I 

| j 

The  next  step  in  the  design  space  selection  was  to  determine  a 
reasonable  range  for  each  variable.  The  ranges  were  chosen  so  as  to 
allow  a second  order  polynomial  in  terms  of  the  Independent  variables 
to  approximate  the  dependent  parameters  (to  be  discussed  in  Chapter  IV). 
The  ranges  in  this  study  were  based  on  previous  aircraft  of  the  same 
type  having  a comparable  mission.  The  range  of  each  variable  was  then 
divided  into  four  equal  segments  so  that  all  four  variables  would  have 
five  possible  values.  The  possible  values  of  the  four  parameters  were 
chosen  as  follows: 


(Airframe  Variable) 

WOS 

70 

80 

90 

100 

110 

(Airframe  Variable) 

AR 

2.5 

3.0 

3.5 

4.0 

4.5 

(Engine  Variable) 

TW 

.5 

.7 

.9 

1.1 

1.3 

(Mission  Variable) 

RN 

150 

200 

250 

300 

350 

where  WOS  has  units  lbs/sq  ft  and  RN  has  units  of  nautical  miles.  There 
are  5^  possible  combinations  that  made  up  the  complete  design  space. 
These  625  design  points  were  input  to  the  CISE  program  which  gave  as 
output  625  TOGW's  and  DISTL's. 

Latin  Square  Design  Selector 

The  computer  time  required  to  run  625  mission  simulations  was 
approximately  600  seconds.  The  need  to  reduce  this  time  and  yet 
maintain  some  degree  of  accuracy  in  the  final  optimal  design  point  led 
to  the  use  of  a statistical  design  selector.  The  latln  square  design 
selector  will  select  out  of  the  entire  design  space  an  appropriate 
subset  which  yields  a good  representation  of  the  complete  space.  The 
representation  is  good  in  the  sense  that  approximating  equations  for 
the  dependent  variables  based  on  the  latin  square  subset  will  predict 
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the  actual  or  trua  valua  to  within  five  percent.  The  staple  latin 
square  technique  uses  a method  (to  be  discussed  later)  which  yields  a 
sequence  of  values  for  each  variable.  The  sequence  Is  formed  by  join* 
ing  together  permutations  of  the  various  values  of  the  variables  in 
such  a way  that  each  variable  le  stepped  through  all  its  values  once 
every  "1M  data  points*  where  "1"  is  the  msaber  of  values  for  each 
variable. 

The  number  of  data  points  selected  by  the  latin  square  technique 
depends  on  the  number  of  independent  variables.  If  there  are  "n"  de- 
sign variables  then  there  will  be  (n*l)*(n*2)/2  terms  in  the  regression 

equation  and  the  regression  routine  will  require  at  least  this  many 

2 

data  points.  The  staple  latin  square  will  generate  (nel)  design  points 
that  have  design  values  of  each  independent  variable  equally  spaced  on 
the  normalised  interval  (-1*1)  (Ref  2:4-5).  Accurate  fits  may  not  be 
obtained  with  this  small  number  of  points  so  the  simple  latin  square 
technique  was  modified  to  generate  more  data  points  (this  was  done  to 
gain  more  accuracy  and  at  the  same  time  tne  time  for  computer  simulations 
of  the  mission  was  increased).  Table  I contains  the  number  of  data 
points  required  by  the  modified  latin  square  for  a given  number  of 
variables. 

The  staple  latin  square  technique  uses  field  algebras*  that  are 
derived  from  modulo  (or  residues  class)  arithmetic  using  the  algebra  of 
integers*  to  establish  "n"  matrices  (where  again  "n"  is  the  number  of 
independent  design  variables)  that  are  used  to  generate  the  design 
points  (Ref  6:57).  Table  II  contains  the  matrices  for  the  four  variable 
case. 
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The  data  spaca  la  obtained  by  finding  the  mid-point  and  an  appro- 
priate step  else  for  each  variable.  The  first  element  of  the  data 
space  is  obtained  by  multiplying  each  variable's  step  sise  by  element 
(1,1)  of  that  variable's  matrix  and  adding  to  the  mid-point.  The 
second  design  point  is  obtained  by  multiplying  each  variable's  step 
site  by  element  (2,1)  and  adding  to  the  mid-point.  This  procedure  Is 
continued  until  the  space  Is  coexists. 

Several  problems  are  possible  when  using  the  simple  latln  square 
technique.  The  most  significant  problem  Is  that  the  selected  set  of 
data  points  may  not  be  Independent  as  shown  In  Fig.  1 (Ref  6t59).  The 
bad  fit  does  not  adequately  represent  the  space,  while  the  good  fit 
does.  The  lack  of  Independence  allows  the  possibility  of  significant 
cross-correlations  building  up  between  the  variables;  however,  these 
cross-correlations  would  be  detected  during  the  surface  fit  procedure. 
These  problems  did  not  appear  in  this  study. 

Inputs  and  outputs  to  the  latln  square  computer  program  are 
described  in  Appendix  A. 


III.  Mission  Simulation 


The  mission  profile  used  in  this  study  is  given  in  Fig.  2 from 
Ref.  5.  This  mission  along  with  the  independent  design  variables  was 
used  as  input  to  a mission  simulation  program  used  by  the  AFFOL  Design 
Branch.  The  CISE  program  was  used  to  calculate  the  TOGW  required  to 
accomplish  the  specified  mission  given  the  set  of  design  variables. 

The  mission  was  based  on  a typical  interdiction  mission  which  an 
advanced  tactical  fighter  might  be  required  to  fulfill.  The  vehicle 
•as  to  perform  the  mission  with  a one  man  crew  and  was  to  carry  5300  lbs 
of  stores.  The  only  requirement  on  the  combat  performance  was  the  fuel 
allowance  of  fifteen  percent  of  the  total  initial  fuel  weight  (Ref  5). 
This  requirement  was  met  by  changing  the  number  of  turns  in  the  combat 
segment  of  the  mission. 

The  CISE  program  was  developed  to  be  used  as  a design  tool  at  the 
"zero  level"  of  design,  i.e.,  even  before  a drawing  of  the  configuration 
is  prepared  (Ref  4:2).  The  program  performs  a weight  oriented  sizing 
process  to  produce  the  basic  physical  characteristics  of  the  vehicle. 

The  procedure  was  set  up  to  use  DO  loops  for  the  four  design  variables. 
An  iterative  process  is  accomplished  for  each  combination  of  variables. 
The  Iteration  assumes  an  initial  TOGW  and  then  calculates  a geometry 
based  on  the  guessed  weight  and  input  variables.  The  input  mission  is 
then  "flown"  and  fuel  needs  are  calculated  to  meet  the  requirements  of 
takeoff  reserves,  combat  penalties,  and  landing  reserves.  The  final 
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process  in  the  iteration  is  to  calculate  the  major  weight  components 
that  include  the  wing,  the  body,  the  tall,  landing  gear,  systems  and 
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1.  Warm  up  and  takeoff 

2.  Cl lob  to  45000  ft 

3.  Cruise  at  .8M  for  200  NM 

4.  Accelerate  to  1.6M 

5.  Climb  to  50000  ft 

6.  Cruise  at  1.6M  for  variable  distances  (150,  200, 

250,  300,  350  NM) 

7.  Combat  (15X  of  initial  fuel  for  combat  penalty) 

8.  Drop  stores 

9.  Cruise  at  1.6M  at  55000  ft  for  same  distance  as  Item  6 

10.  Decelerate  to  ,8M  and  cruise  back  at  45000  ft 

11.  Loiter  at  sea  level  for  20  minutes  and  land 

Fig.  2.  Mission  Profile 
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equipment,  and  payload.  These  group  weights  are  sunned  to  yield  a TOGW 
at  the  end  of  the  Iteration  cycle.  When  this  weight  Is  within  a speci- 
fied tolerance  of  the  Initial  guess,  then  the  process  Is  terminated  and 
another  design  point  is  Initiated.  If  the  weight  Is  not  within  the 
tolerance,  then  it  Is  used  as  the  Initial  guess  and  the  process  is 
repeated. 

The  procedure  that  would  normally  be  followed  In  the  design  process 
using  surface  fit  approximations  would  be  to  Input  the  latln  square  de- 
sign space  of  independent  variables  (in  this  particular  study  the  45 
design  points  of  Table  III)  to  the  mission  simulator  which  would  output 
the  dependent  variables  (In  this  particular  study  the  TOGW's  and  DISTL's 
of  Table  IV).  These  steps  would  be  followed  by  a surface  fit  procedure 
and  an  optimisation  routine  (both  to  be  discussed  In  later  chapters). 
However,  this  normal  procedure  was  deviated  from  slightly  in  this  study. 

It  was  desired  to  have  a basis  for  judging  the  "goodness  and  accuracy" 
of  the  surface  fit  approximations.  In  order  to  accomplish  this  goal, 
the  entire  design  space  consisting  of  all  625  design  points  were  Input 

s 

to  the  CISE  program  which  output  625  TOGW's  and  DISTL's.  Those  points 
which  had  been  selected  by  the  latln  square  technique  (Table  III)  were 
included  in  the  625  design  points;  hence,  the  45  output  functions 
(TOGW's  and  DISTL's)  for  the  latln  square  design  space  were  already 
available  from  having  run  the  simulation  for  the  complete  design  space. 

There  would  have  been  a slight  difference  In  the  values  of  the 
output  functions  for  the  45  latln  square  design  points  had  they  been 
the  only  data  points  Input  to  CISE  as  Is  the  case  In  the  normal  design 
procedure  described  above.  This  is  due  to  the  weight  tolerances  and 
Initial  guess  in  the  Iteration  process.  These  differences  were  negligible. 


The  TOGW's  for  the  selected  design  space  are  given  In  Table  III 


In  the  results  section. 

Input  and  output  requirements  for  the  CISE  program  are  given  In 
Appendix  B. 


IV.  Surface  Fit  Procedure 


The  next  step  In  the  process  was  to  perform  a regression  analysis 
on  the  performance  function  data.  A second  order  polynomial  was  to  be 
generated  relating  the  TOGW  to  the  four  design  variables,  l.e.. 


n n n 

TOGW  ■ A-  ♦ E A1X1  ♦ E E Aii  X*  X« 

1-1  j-1  1-1  J J 


where  the  A's  are  the  coefficients  and  the  X's  are  the  four  design 
variables.  The  design  variables  were  Input  In  a normalized  form  so 
that  all  values  would  lie  on  the  Interval  (-1,1).  This  was  done  by 


using  the  following  transformation: 


x - % (Xmax  ♦ Xm1n) 

* (XMAX  " XMIN* 


This  normalization  was  necessary  to  eliminate  round-off  errors  which 
would  have  occurred  due  to  the  difference  In  magnitude  of  the  variables. 
The  normalized  variables  were  also  used  in  the  optimization  routine  and 
transformation  to  real  range  values  occurred  when  an  optimal  design 
point  had  been  found. 

The  quadratic  approximation  was  chosen  for  two  basic  reasons.  The 
first  reason  was  that  the  quadratic  polynomial  is  simple  and  easy  to 
work  with  and  makes  possible  the  easy  and  economical  calculation  of 
partial  derivatives  used  in  the  optimizer.  Adequate  representation 
can  be  obtained  with  second  order  surfaces  thus  avoiding  higher  order 
approximations  (Ref  7:595).  The  second  reason  is  more  mathematically 
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rational.  Basically,  the  second  order  assumption  Is  equivalent  to 
assuming  that  the  performance  function  can  be  adequately  approximated 
by  the  first  few  terms  of  Its  Taylor  series  expansion  about  some  nominal 
design  point  (x^,  x2,  Xj,  x^).  The  basic  requirement  Is  that  the  range 
of  values  for  the  design  variables  be  kept  reasonably  small.  If  this 
Is  the  case,  then  the  quadratic  polynomial  derived  from  the  regression 
analysis  will  provide  adequate  approximations  within  these  ranges 
(Ref  2:4). 

The  number  of  unknown  regression  coefficients  (K)  In  the  quadratic 
polynomial  of  Eq  (1)  Is  (n+l)*(n*2)/2  where  "n"  Is  the  number  of  Inde- 
pendent variables.  The  number  of  design  data  points  must  be  greater 
than  or  equal  to  the  number  of  unknown  coefficients  so  that  an  over- 
determined  system  of  linear  equations  In  the  form  of  Eq  (1)  In  terms 
of  the  unknown  coefficients  Aq,  Aj,  and  Ajj  can  be  solved  by  the  method 
of  least  squares.  The  least  squares  method  selects  the  terms  so  as  to 
minimise  the  error  (or  residual)  sun- of -square,  SSE, 

N 2 

SSE  - E €„  (3) 

l-l  * 

where  €R  Is  the  difference  between  the  true  value  of  the  performance 
function  and  the  approximated  value  and  where  N Is  the  number  of  data 
points  Input  to  the  regression  analysis. 

The  goodness-of-f It  of  the  regression  surface  Is  tested  statis- 
tically by  variance  analyals.  The  regression  program  used  In  this 
study  employed  four  tests  to  determine  the  goodness-of-f It.  These  were 
the  standard  F-statlstlc  for  regression,  the  multiple  correlation  co- 
efficient squared,  the  significance  ratio,  and  the  atandard  error. 
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The  F-statlstlc  is  defined  as  the  ratio  of  the  regression  mean 


1 


square, 

N _ 2 

Z (TOGWjj  - TOGW) 

MSR  - k"1  (4) 

N 

to  the  error  mean  square, 

N 2 

Z (TOGWk  - TOGW.) 

MSE  - ■ (5) 

(H  - K) 

where  TOGW  is  the  mean  of  the  true  TOGW’s  in  the  N data  points,  TOGWk 
is  the  TOGW  predicted  by  the  polynomial  at  data  point  k,  and  K is  the 
number  of  coefficients  in  Eq  (1).  The  F value  calculated  should 
exceed  the  F value  found  in  the  standard  tables  for  K and  N-K  degrees 
of  freedom  at  the  95%  confidence  level  to  ensure  a good  fit  (Ref  6:595). 
If  this  test  is  met,  then  a significant  regression  is  assumed. 

The  multiple  correlation  coefficient  squared  (MCC2  or  R2)  is 
defined  as 

N _ 2 

Z (TOGWk  - TOGW) 

MCC2  - — (6) 

Z (TOGWj.  • TOGW)2 
to-1 

where  the  terms  are  as  defined  for  Eqs  (4)  and  (5).  This  quantity 


The  standard  error  Is 


STERR 


(7) 


and  is  a measure  of  goodness-of-f It.  The  smaller  the  value  of  the 
standard  error  the  better  the  approximation. 

Information  on  the  significance  ratio  used  in  the  regression  pro- 
gram was  limited  to  the  fact  that  it  should  be  as  large  as  possible  in 
order  to  ensure  a good  fit.  The  significance  ratio  was  used  to 
determine  which  step  in  the  regression  analysis  would  be  used  to  select 
the  coefficients  for  the  quadratic  approximation.  A forward  and  back- 
ward regression  analysis  was  performed,  i.e.,  variables  were  entered 
and  removed  whenever  a better  fit  could  be  obtained.  The  second  order 

and  cross  product  terms  were  considered  as  separate  variables  in  this 

2 

part  of  the  problem;  that  is,  (WOS  ) could  be  entered  at  a step  even 
though  WOS  was  not  in  the  equation.  The  step  chosen  to  be  the  most 
representative  was  based  on  maximizing 


FACT 


I SIG  RATIO 
J STERR 


MCC2 


(8) 


Input  and  output  requirements 


for  the  regressor  are  found  in 


Appendix  C 


Surface  Fit  Optimization 

The  basic  optimization  problem  considered  was  to  minimize  or 
maximize  some  performance  function  subject  to  constraints  on  other 
performance  functions.  A maximization  problem  can  be  converted  to  a 
minimization  problem  simply  by  taking  a negative  of  the  function  to  be 
maximized  and  then  performing  a minimization  (l.e.,  RANGE  Is  maximized 
by  minimizing  -RANGE).  Similarly,  inequality  constraints  can  all  be 
regarded  as  less  than  or  equal  to  constraints.  This  is  accomplished 
by  taking  the  negative  of  a function  that  is  to  be  greater  than  some 
value  and  making  it  less  than  the  negative  of  this  value  (i.e., 

RANGE  > 150  would  be  rewritten  -RANGE  < -150  where  -RANGE  would  be 
defined  as  a new  function).  Thus,  using  these  conventions,  the 
optimum  design  problem  can  be  rewritten  as 

Minimize  fj  (x) 

Subject  to  fj  (x)  < (l-2,...,m)  (9) 

where  x is  the  vector  of  independent  variables,  f's  are  the  perform- 
ance functions,  and  the  c's  are  the  fixed  numerical  values  defining 
the  upper  limits  on  the  constraint  functions.  The  Independent 
variables  are  also  constrained  to  lie  within  a region  of  interest, 
within  which  an  optimal  design  is  supposed  to  lie.  This  region  is  the 
area  determined  in  the  selection  of  the  design  space  and  is  the  region 
in  which  the  quadratic  approximation  is  valid. 

The  penalty  function  method  has  been  shown  to  work  on  the  minimi- 
zation problem  of  Eq  (9)  (Ref  6:29-32).  This  method  transforms  the 
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inequalities  out  of  the  problem  by  establishing  a penalized  cost 
function  of  the  form 

m , 

F(x)  - f j(x)  ♦ PK  E^  (CVj r (10) 

where  CV  represents  the  violation  of  the  constraint  Inequalities;  if 
the  constraint  is  satisfied  then  CV  » 0.  PK  Is  a weighting  factor  in 
Eq  (10)  which  modulates  the  severity  of  violating  the  constraints. 

This  weighting  factor  can  be  considered  the  tolerance  level  for  con- 
straint violations;  hence,  if  the  acceptable  deviation  from  the  con- 
straints was  one  percent  then  PK  would  be  100.  It  was  assumed  in  this 
study  that  all  constraints  required  this  tolerance  and  PK  was  taken  to 
be  a constant  equal  to  100. 

The  penalty  function  method  has  the  advantage  that  the  non-linear 
constraints  are  removed  from  the  optimization  problem  by  a simple 
transformation  approach.  It  also  has  an  advantage  in  that  it  is  not 
necessary  to  locate  a feasible  starting  point  for  the  optimization 
search.  The  penalty  function  algorithm  can  always  be  started  from  a 
ncnfeaslble  point  (Ref  6:31). 

There  are,  however,  two  major  problems  in  the  penalty  function 
method.  The  first  problem  is  that  for  large  values  of  the  weighting 
factor  PK  in  Eq  (10),  the  penalty  cost  function  can  develop  a long 
narrow  valley  containing  the  minimum  point  (Ref  1:40-41).  Gradient 
procedures  for  finding  the  minimum  tend  to  go  back  and  forth  between 
the  sides  of  the  valley  instead  of  going  down  the  long  direction  of 
the  valley.  Another  problem  is  that  the  generation  of  the  penalty 
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function  My  also  create  new  pseudo  minlM  that  are  not  present  In  the 
original  problem.  Ref  2 (p.  31)  states  "These  difficulties  are 
numerical  in  nature,  corresponding  to  large  values  of  PK  and/or  a 
sharply  decreasing  cost  function.  These  difficulties  should  not  occur 
with  the  approximating  functions  which  are  used." 

The  form  of  the  quadratic  functions  supplied  by  the  regression 
analysis  was 

f|(jc)  - Ajo  ♦ Ajj  Qc)  ♦ xT  Aj2  x (1-1 m)  (11) 


where  A^  Is  a constant  term,  A^  is  an  "n”  vector,  and  A^  Is  an  n x n 
symmetric  Mtrlx.  The  form  of  this  equation  suggests  the  use  of  a 
gradient  type  search  procedure  since  the  gradient  of  the  cost  function 
would  be  readily  available.  Indeed,  the  desire  to  use  a gradient 
method  was  one  of  the  reasons  for  attempting  the  surface  fit.  The 
gradient  of  the  penalty  cost  function  Is 


g (x ) 


m 

A11  * a12  & ♦ 2*k 


(cv j ) (Aj  i ♦ Aj2  x) 


(12) 


A conjugate  gradient  technique  was  used  in  this  study.  The 
ordinary  gradient  has  the  drawback  of  slow  convergence  and  methods 
have  been  developed  to  reduce  the  number  of  iterations  required  for 
convergence  (Ref  9:10).  These  methods  use  a modified  search  direction 
to  speed  up  the  convergence  rate. 

The  minimising  algorithm  used  was  as  follows:  (a)  an  evaluation 

of  the  gradient  (gj)  at  a given  point  xj  was  accomplished,  (b)  the 
search  direction  pj  was  computed  according  to 


! 


20 


Pi  - Vl, 


where  "i"  denotes  the  present  point  and  yas  given  by 


Hi  - Himl  ♦ 


<A  * i-1  - Hi-l  A (A  xlml  - Ht!t  A g^)1 

<A  *1-1  " H,Ii  A *,.i>  Ait.! 


where  1-1  denotes  the  previous  point,  Hq  was  defined  to  be  the  unity 
matrix,  - gj  - g^,  and  Xj.i  • Xj  - x^j*  * one  dimensional 

search  was  performed  on  the  function  FQc^j)  " F(*i  + AjCj)  where 

^1  " "ti  £i  (*5) 


and  the  search  was  to  minimize  F(xj  - tj  £j)  with  respect  to  tj  (this 
one  dimensional  search  will  be  discussed  later),  (d)  having  found  the 
t}  that  yielded  a minimum,  the  displacement  x^  was  computed  according 
to  Cq  (IS),  (e)  the  next  point,  x^j  was  obtained  by 

&I+I  " *1  ♦ A*j  <16^ 

The  value  of  will  always  be  less  than  the  value  of  F(Xj) 

(Ref  8:271),  These  steps  (a-e)  were  repeated  until  the  stopping 


condition 


Si*i  £ 10’6 


was  met.  This  stopping  condition  may  correspond  to  either  a minimal 
point  or  a nonmlnimal  stationary  point  (Ref  8:275),  For  this  reason, 
small  perturbations  about  this  point  were  Investigated. 


i 


t 


The  entire  algorithm  may  need  to  be  restarted  by  setting  - HQ. 
This  Is  required  when  the  Inequality 


Wi  pi 

Is  satisfied.  If  the  algorithm  Is  not  restarted  on  this  condition 
then  the  starting  point  for  the  one-dimensional  search  will  automatic- 
ally  satisfy  the  search  stopping  condition  of  Eq  (22)  so  that  tj  - 0. 

The  linear  search  method  was  that  proposed  by  Davldon  and  has  been 
used  by  Fletcher  and  Powell  (Ref  3).  The  method  was  (*)  to  estimate 
the  order  of  magnitude  of  the  stepsise  tj,  (b)  establish  bounds  on  t 
between  which  the  minimus  point  lies,  and  (c)  to  perform  a cubic 
interpolation  to  find  tj. 

The  determination  of  the  order  of  magnitude  of  tj  was  based  on  an 
available  estimate  (est)  of  the  value  of  FQc)  at  the  unconstrained 
minimum.  This  estimate  was  available  in  the  problem  at  hand  and  would 
be  expected  to  be  around  -1  when  normalised  variables  were  being  used. 
The  step  sise  was  determined  by 


10 


•6 


(18) 


or 


h - k 


2(F(Xi)  - est) 
T 

Pi  *1 


if  0 < k < 


h 


4* 


otherwise. 


The  bounds  on  tj  were  established  by  examining 


y<t)  - FQcj  - t pj) 
y'(t)  - -£i  4 <&i  • t £|> 


(20) 
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at  points  where  t ■ 0,  h,  2h,  4h,  ... , a,  b where  b is  the  first  point 
at  which  either  y’  la  non-negative  or  y has  not  decreased.  It  then 
follows  that  tj  Is  bounded  In  the  Interval  (a,b). 

Given  these  Halts  on  tj  a cubic  interpolation  estimating  tj  may 
be  given  by 


t_  - b - ( y*  (b)  ♦ w - t v (b  - a) 

• 1 y’ (b)  - y* (a)  ♦ 2w  ' * 


* - 3 < y-W.  " VW  ) ♦ y.(a)  ♦ 


b • a 


y'  (b) 


tr  - y* (a)  y'(b) 


If  neither  y(a)  nor  y(b)  Is  less  than  y(te)  then  te  Is  accepted  as  the 
estimate  of  tj  and  the  stopping  condition  is  tested  at  te.  If  the 
condition  Is  met,  the  one  dimensional  search  is  terminated.  If  the 
condition  is  not  met,  the  cubic  Interpolation  Is  repeated  over  the 
sub-interval  (a,te)  if  yf(tft)  is  positive  and  the  cubic  interpolation 
is  repeated  over  the  sub-interval  (t#,b)  if  y'(t4)  Is  negative.  The 
stopping  condition  Is  determined  by  precision  requirements  and  was 


Slel  Ei  * l0'6 


where  the  constraint  limit  must  be  the  same  as  in  Eq  (18). 

Two  different  minimisation  problems  were  investigated.  The  first 
was  to  minimise  the  TO GW  subject  only  to  the  constraints  on  the  inde- 
pendent design  variables.  In  solving  this  problem,  the  design 
variables  were  entered  into  the  penalty  cost  function  by  'creating 
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Cals*  performance  functions  consisting  only  of  the  variables  and  their 
limits.  The  second  problem  was  to  minimize  the  TOGW  subject  to  a 
constraint  on  the  landing  distance  which  was  a performance  parameter 
available  from  CISE. 

Input  and  output  requirements  for  the  optimization  routine  are 
described  in  Appendix  0. 


VI.  Results 

The  design  space  selected  by  the  Latin  Square  technique  is  given 
in  Table  III.  Each  design  point  was  input  to  the  mission  simulator  and 
the  performance  values  (TOGW  and  OISTL)  for  each  point  are  given  in 
Table  IV.  These  45  points  were  then  input  to  the  regression  analysis. 
The  second  order  polynomial  generated  for  each  performance  function 
was  (equations  are  in  terms  of  normalized  values) 

TOGW  - -.31638  ♦ .06307  (WOS)  - .04207  (AR)  ♦ .79413  (TW) 

♦ .18533  (RN)  ♦ .03678  (WOS2)  - .04107  (WOS)(AR) 

- .04316  (WOS) (TW)  - .03023  (AR)(TW)  - .31439  (TW2) 

♦ .07082  (TW) (RN) 

(23) 

OISTL  - .01156  ♦ .87538  (WOS)  - .10542  (AR)  - .05551  (TW) 

- .08192  (RN)  - ,02251  (WOS2)  - .01018  (W0S)(AR) 

- .01253  (W0S)(TW)  - .02614  (WOS)(RN)  - .00507  (AR2) 

♦ .00966  (AR) (TW)  ♦ .00488  (AR) (RN)  - .02229  (TW2) 

♦ .00287  (TW) (RN) 

The  statistical  analysis  for  each  equation  is  given  in  Table  V. 

The  value  of  the  P-statistic  from  the  standard  tables  for  a numerator 
of  45  degrees  of  freedom  and  a denominator  of  15  degrees  of  freedom  is 
2.19.  The  TOGW  equation  does  meet  this  requirement  so  it  has  a 95% 
confidence  level  that  it  represents  the  true  data  space.  The  other 

1 

(values  for  the  TOGW  equation  also  Indicate  a good  fit.  The  P-value 

for  the  OISTL  equation  was  lower  than  the  required  value;  however,  the 
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Table  III 

Selected  Latin  Square  Design  Space 


I 


Case  No.  WOS 


1 90 

2 100 

3 110 

4 70 

5 80 

6 110 

7 80 

8 100 

9 70 

10  100 

11  110 

12  70 

13  80 

14  90 

15  70 

16  90 

17  110 

18  80 

19  110 

20  70 

21  80 

22  90 

23  100 

24  80 

25  100 

26  70 

27  90 

28  70 

29  80 

30  90 

31  100 

32  110 

33  90 

34  110 

35  80 

36  100 

37  80 

38  90 

39  100 

10 
70 
00 

>3  70 

4 90 

>5  110 


AR 

TW 

RN 

3.5 

.9 

250 

4.5 

.5 

200 

3.0 

1.1 

150 

4.0 

.7 

350 

2.5 

1.3 

300 

4.0 

.7 

150 

4.5 

.5 

300 

2.5 

1.3 

200 

3.0 

1.1 

350 

4.0 

1.1 

300 

2.5 

.7 

250 

3.5 

1.3 

200 

4.5 

.9 

150 

3.0 

.5 

350 

4.5 

.9 

200 

2.5 

.7 

350 

3.0 

.5 

250 

3.5 

1.3 

150 

4.5 

1.3 

350 

3.0 

.9 

300 

4.0 

.5 

250 

2.5 

1.1 

200 

3.5 

.7 

ISO 

2.5 

1.1 

250 

3.0 

.9 

150 

3.5 

.7 

300 

4.0 

.5 

200 

2.5 

.5 

150 

3.5 

1.1 

350 

4.5 

.7 

300 

3.0 

1.3 

250 

4.0 

.9 

200 

3.0 

1.3 

300 

3.5 

1.1 

200 

4.0 

.9 

350 

4.5 

.7 

250 

3.0 

.7 

200 

4.0 

1.3 

150 

2.5 

.9 

350 

3.5 

.5 

300 

4.5 

1.1 

250 

3.5 

.5 

350 

4.0 

1.3 

250 

4.5 

1.1 

150 

2.5 

.9 

300 

■ 
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Table  IV 

Mission  Simulation  Output 


Case  No. 

TOGW 

0ISTL 

1 

35407 

6230 

2 

31751 

6831 

3 

38158 

7734 

4 

33579 

4721 

5 

45299 

5514 

6 

32253 

7655 

7 

31567 

5380 

8 

43695 

6981 

9 

40411 

4812 

10 

39904 

6671 

11 

33844 

7745 

12 

42660 

4812 

13 

33956 

5488 

14 

32604 

6246 

IS 

34885 

4746 

16 

33936 

6307 

17 

33302 

7690 

18 

41480 

5560 

19 

45256 

7091 

20 

36169 

4879 

21 

31117 

5522 

22 

38776 

6393 

23 

31597 

7096 

24 

39253 

5655 

25 

34302 

7144 

26 

32873 

4852 

27 

31240 

6258 

28 

30071 

5136 

29 

40393 

5409 

30 

32821 

6033 

31 

44253 

6837 

32 

34889 

7530 

33 

44842 

6129 

34 

38751 

7550 

35 

36755 

5354 

36 

32248 

6763 

37 

31852 

5727 

38 

41383 

6152 

39 

37548 

6916 

40 

33546 

7513 

41 

38687 

4691 

42 

33215 

6805 

43 

43414 

4703 

44 

37176 

6124 

45 

37297 

7625 
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other  statistical  values  are  extremely  good  and  the  maximum  difference 
that  the  equation  predicts  from  the  true  data  Is  only  six  percent. 


Table  V 

Statistical  Analysis  of  Selected  Equations 


P-value 

Significance 

Ratio 

Multiple 

Correlation 

Standard 

Error 

TOGW 

5.3352 

922.4962 

.99616 

.04042 

OISTL 

.8058 

25053.265 

.99995 

.00758 

One  problem  was  noted  In  the  TOGW  equation.  The  AR  appeared  only 
as  a linear  term  Indicating  that  a proper  minimum  would  not  exist. 
Also,  the  sufficiency  condition  for  a minimum  was  not  met  in  that  the 
matrix  of  second  partial  derivatives  was  not  positive,  semi-definite. 

To  analyse  this  problem,  graphs  were  made  of  the  TOGW  versus  the 
AR  for  several  values  of  WOS  at  constant  RN  and  TW.  These  graphs  are 
presented  In  Pigs.  3-8.  It  was  concluded  from  these  graphs  that  the 
calculated  curves  do  represent  the  actual  data  trends  and  yield  values 
of  TOGW  that  are  within  two  percent  of  the  actual  data. 

Three  dimensional  plots  were  made  of  TOGW  versus  AR  and  WOS  at 
constant  RN  and  TW.  Once  again  the  calculated  plots  exhibit  the  same 
shaped  (although  they  are  smoother  as  would  be  expected)  as  the  true 
data  and  have  a maximum  difference  of  two  percent.  These  surfaces  are 
presented  In  Pigs.  9-14. 

Two  minimisations  were  performed.  The  first  problem  was  to  find 
the  minimum  TOGW  within  the  design  space.  The  second  problem  was  to 
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find  the  minimum  TOGW  subject  to  the  constraint  that  the  OISTL  be  less 
than  5000  feet  and  still  be  within  the  design  space.  These  minlmums 
were  also  found  by  manually  analysing  the  625  data  points  that  repre- 
tented  the  complete  data  space.  The  minimum  point  for  each  range  was 
found  and  the  results  are  presented  in  Tables  VI  and  VII. 

Table  VI  indicates  that  the  surface  fit  approximations  yield 
minimum  values  of  TOGW  that  are  within  two  percent  of  the  actual  mini* 
mum  values.  The  design  parameters  are  identical  with  the  exception  of 
the  AR.  The  graphs  of  Figs.  3-8  indicate  that  the  TOGW  is  not  signifi- 
cantly affected  by  changes  In  A&  so  the  problem  of  different  AR's  is 
acceptable.  Table  VII  shows  similar  results  for  the  case  when  OISTL 
was  constrained  to  be  less  than  or  equal  to  5000  feet. 

These  results  also  indicate  that  the  minimum  design  point  is  not 
Inside  the  range  of  the  design  variables.  In  the  unconstrained  problem, 
the  design  point  was  at  the  lower  limit  of  each  of  the  variables.  This 
does  present  a problem  in  that  the  minimum  was  assumed  to  lie  within 
the  design  region.  However,  the  approximation  method  did  predict  the 
same  value  as  was  found  by  looking  at  the  complete  space. 


TOGW  *10 

280-0  285-0  290-0  295-0  300-0  30S-0  310-0  315-0  320-0  325-0  330-0  335-0  340-0 


TOGW  VS  flR 

RN  - 150,  TW  - .5 
CALCULATED  VALUES 


LEGEND 
□ -WOS  - 70 
o-WOS  - 80 
a-WOS  “ 90 
+ -WOS  - 100 
x - WOS  - 110 


Ftg.  4.  TOGW  vs.  AR  (RN  - 150,  TW  - .5)  - Calculated  Values 
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TOGW  VS  RR 


34*0 


TOGW  VS  WOS  - RR 

RN  - 150,  tw  - .5 
ACTUAL  DATA 


Pig.  9.  TOGW  v*.  WOS  - AR  (RN  - 150,  TW  - .5)  Actual  Data 
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TOGW  VS  WOS 


RR 


RN  - 150,  TW  - .5 
CALCULATED  VALUES 


rig.  10. 


TOO)  »•.  WOS  - AK  (RH  . 150,  TU  . .5)  Calculated  Values 


Fig.  11.  TOGV  v*.  W03  « AR  (RN  » 250,  TW  - .9)  Actual  Oita 


•if 


44*0  V5*0  48*0  47*0  48*0 


TOGW  VS  NOS  - OR 


RN  - 350,  TW  - 1.3 
CALCULATED  VALUES 


Table  VI 

Results  of  the  Minimisation  - DISTL  Not  Constrained 


Case 

RN 

WOS 

AR 

TV 

TOGW 

A 

150 

70 

3.5 

.5 

29864 

B 

150 

70 

2.5 

.5 

30000 

A 

200 

70 

3.0 

.5 

30295 

B 

200 

70 

2.5 

.5 

30436 

A 

250 

70 

3.0 

.5 

30774 

B 

250 

70 

2.3 

.5 

30872 

A 

300 

70 

3.0 

.5 

31282 

B 

300 

70 

2.5 

.5 

31308 

A 

350 

70 

3.0 

.5 

31818 

B 

350 

70 

2.5 

.5 

31743 

where  Case  A represents  the  complete  data  space  analysis;  Case  B repre- 
sents the  surface  fit  approximation  analysis. 


Table  VII 

Results  of  the  Minimisation  • DISTL  Constrained 


Case 

RN 

WOS 

AR 

TV 

TOGW 

DISTL 

A 

150 

70 

3.5 

.5 

29864 

5007* 

B 

150 

70 

3.5 

.5 

30223 

4986 

A 

200 

70 

3.5 

.5 

30301 

4947 

B 

200 

70 

3.5 

.5 

30659 

4941 

A 

250 

70 

3.0 

.5 

30773 

4965 

B 

250 

70 

3.0 

.5 

30983 

4974 

A 

300 

70 

3.0 

.5 

31282 

4905 

B 

300 

70 

2.5 

.5 

31308 

5000 

A 

350 

70 

3.0 

.5 

31816 

4844 

B 

350 

70 

2.5 

.5 

31744 

4951 

where  Case  A represents  the  complete  data  space  analysis; 
sents  the  surface  fit  approximation  analysis. 

♦This  number  Is  within  the  constraint  tolerance. 

Case  B rep re- 

VH.  Conclusions  and  Recommendations 


Conclusions 

1.  The  surface  fit  approximations  can  be  used  in  the  design 
process  to  predict  designs  that  agree  with  more  detailed  and  time  con- 
suming processes. 

2.  The  fact  that  the  independent  variables  are  constrained  to 
lie  in  the  selected  region  requires  that  care  be  taken  during  the 
definition  of  the  ranges  of  these  variables.  The  ranges  must  be  chosen 
small  enough  so  that  the  quadratic  approximations  are  valid;  however, 
they  must  be  large  enough  to  allow  the  minimum  point  to  lie  within  the 
design  space.  This  will  prevent  a minimum  design  point  which  satisfies 
the  performance  function  constraints  from  not  being  found. 

3.  The  simple  latln  square  technique  as  modified  will  yield  a 
representative  design  space.  This  reduces  significantly  the  number  of 
mission  simulations  that  must  be  performed;  hence,  computer  time  and 
engineering  analysis  are  reduced. 

Recommendations 

1.  Further  studies  should  be  made  to  demonstrate  the  capabilities 
with  an  expanded  number  of  variables  and  performance  functions.  These 
studies  should  be  made  in  close  coordination  with  the  AFFDL  Design 
Branch  using  a design  problem  that  they  are  considering. 

2.  Other  design  space  selection  techniques  such  as  the  orthogonal 
latin  square  and  the  D- optimal  should  be  considered.  These  techniques 
may  further  reduce  the  number  of  mission  simulations. 
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3.  It  la  appropriate  to  Incorporate  the  surface  fit  approximations 
In  the  design  procedures  at  APFOL.  These  techniques  are  currently  being 
used  by  some  aerospace  companies  and  In  particular  at  Boeing,  enabling 
them  to  save  many  computer  and  engineering  hours. 
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Appendix  A 


Latin  Square  Computer  Program  Input/Output 


Card  No,  Col.  No. 
1 12-13 

2-N  1-6 

11-20 

21-30 


Description 

Number  of  independent  design  variables 
(Format  12) 

Alphanumeric  label  for  the  variable 
(Format  A6) 

Lower  limit  for  the  variable 
(Format  CIO. 4) 

Upper  limit  for  the  variable 
(Format  C10.4) 


”N"  is  the  number  of  independent  design  variables  plus  one  which 
allows  the  range  of  each  variable  to  be  input. 

Input  is  assigned  to  Tape  5. 


Output 

The  format  of  the  output  is  (18,  10F12.4/14X, 7F12.4)  which  yields 
the  data  point  number  and  the  values  of  the  design  variables  for  that 
data  point. 

Output  is  assigned  to  Tape  6. 
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Appendix  B 


i.: 


! 


Computerized  Initial  Sizing  Program  (CISE)  Input/Output 

All  Input  is  in  the  standard  Fortran  Namelist  format  with  the  first 
character  in  each  record  appearing  in  the  second  column.  The  title  of 
the  design  under  study  appears  on  the  first  data  record  and  is  in 
Hollerith  format.  The  input  is  divided  into  five  groups  using  the 
following  names: 

DESIGN,  MISSION,  WEIGHTS,  GEOM,  PROP 
The  items  that  can  be  input  by  each  Namelist  are: 

$DESIGN 

LF  - Load  Factor 

VMAX  * Maximum  Equivalent  Airspeed 
AhMAX  - Maximum  Mach  No. 

ALTX  - Altitude  for  AMMAX 

PS  - Energy  Level  Required  at  Flight  Design  Gross  Wt. 

GPS  - G Level  for  PS 

ALTPS  - Altitude  for  PS 

AMNPS  - Mach  No.  for  PS 

NCREW  -No.  of  Crew  Members 

NTANK  - No.  of  300  Gal  External  Tanks 

NSTORX  - No.  of  External  Stores 

NPYL  - No.  of  Pylons  to  Carry  External  Stores  and  Tanks 
NS TOR I - No.  of  Internal  Stores 

CiMAX  - Maximum  Lift  Coefficient  for  Takeoff  and  Landing 
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NCOK  • Cd^  Calculation  Cue;  When  - 0,  CDOSF  » Cf  ; When  - 1, 

Cd  la  calculated  by  Program  and  CDOSF  Is  the  Modification 
o 

Factor 

CDOSF  - Equivalent  Skin  Friction  Coefficient  or  Ch  Modification 

o 

Factor  (depending  on  NCDX  value  input) 

CDS 1X1  • Subsonic  Store  Drag  Modification  Factor 
CDSTX2  - Supersonic  Store  Drag  Factor 
ALTTOL  - Altitude  for  Takeoff  and  landing 

DTEMP  - Takeoff  and  Landing  F from  Standard  Day  at  ALTTOL 
IFROP  ■ 1 for  Turbofan  (P&W  401) 

■ 2 for  Turbojet  (GE  J79) 

■ 3 for  Reciprocating  (Lycoming  IGO-540) 

AB  - 1 for  afterburner 

AB  4 1 for  no  afterburner 

TOWE  - Engine  Thrust  to  Weight  Ratio 

FDGWF  - Ratio  of  Total  Fuel  on  Board  to  Define  Flight  Design  Gross 

Weight  (FDGW)  where,  in  the  Program  FDGW  - TOGW  - (1.0  - FDGWF) 
* WFUEL  (Initialised  at  .60) 

LDGWFS  ■ Same  as  FDGWF,  except  for  Structural  Design  Landing 
Gross  Weight  (LDGW)  where,  in  the  Program 
LDGW  - TPGW  - (1.0  - LDGWFS)  * WFUEL  - WSTOR 
(Initialised  at  .90) 

LDGWF  » Same  as  LDGWFS,  except  for  Landing  Distance  Calculations, 

Using  Design  Landing  Distance  Gross  Weight  (LDGWLD)  where, 
in  the  Program  LDGWLD  - TOGW  - (1.0  - LDGWF)  * WFUEL  • WSTOR 
(Initialised  at  .70) 
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IPRINT  • 1 for  Error  Checkout  Messages  to  be  Printed 


■ 0 for  Error  Checkout  Messages  not  to  be  Printed 
TOTHF  - Thrust  Modification  Factor  for  Non-standard  Day  (Altitude  and 

Temperature  Adjustment)  (Initialised  at  1.0) 

CDP  - Total  CD  Modification  Factor*  May  be  Used  to  Adjust  L/D  by 
Modifying  CD 

COSF  m Factor  to  Adjust  Supersonic  Drag  Contribution  to  CDq 
^MISSION 

N m Total  Number  of  Mission  Legs 
R(I)  - 0 for  Warmup  Fuel  Allowance 
m -l  to  Drop  Stores 

- -2  to  Turn  at  Fixed  Gee 
• -3  to  Accelerate 

■ -4  for  Energy  Altitude  Combat  Fuel  Allowance 

- -5  for  Loiter  Performance 

- -6  for  Loiter  at  Fixed  Altitude  and  Mach  No. 

■ -7  to  Turn  at  Maximum  Possible  Gee 
> 0 for  Climb  Performance 

- 1.E6  to  Climb  on  Intermediate  (Military  Power) 

- 1.1E6  to  Climb  with  Afterburner  with  Distance  Credit 
GEE(I)  - G Level  for  Turns  (no  lg  turns)*  When  RANGE(I)  - -7 

■ Time  (min)  on  Afterburner  Power*  When  RANGE(I)  - 1.1E6 

■ No.  of  External  Stores  to  be  Dropped*  When  RANGE ( I ) - -1 
NTURNS(I)  - No.  of  Turns*  When  RANGE(I)  - -7 

- Mach*  When  RANGE(l)  - -3 

• Energy  Altitude  for  Combat  Fuel  Allowance  (ft)*  When 
RANGE(I)  - -4 
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TIME (I)  • Loiter  Tine  (min).  When  RANG£(1)  - >6,  or  -5 

■ No.  of  Internet  Stores  to  be  Dropped,  When  RANGE(I)  * *1 
I TANK (I)  - No.  of  External  Fuel  Tanks  to  be  Dropped,  When  RANGE(I)  - -1 
AMACH(I)  - Leg  Mach  No.,  Except  When  RANGE(I)  - -1 

■ Weight  of  Cargo  to  be  Dropped,  When  RANGE(I)  - >1 
ALT(I)  - Leg  Altitude 

RG(I)  - Distance  Covered  in  a Particular  Segment,  When  RANGE(I)  ■ 1.1E6 
$ WEIGHTS 

WAVUN  - Weight  of  Avionics  Equipment 

WMA  - Weight  of  AMMO,  Guns,  etc.  That  Does  Not  Require  Installation 
Weight  In  Addition 
WSTORX  - Weight  of  External  Stores 
WSTORI  - Weight  of  Internal  Stores 
WINGF  - Wing  Weight  Modification  Factor 
TAILF  - Tail  Weight  Modification  Factor 
BODYF  m Body  Weight  Modification  Factor 
GEARF  m Gear  Weight  Modification  Factor 

$GEOM 

NWOS  m No.  of  Wing  Loadings  to  be  Cycled,  Maximum  of  S 
WOS(I)  - Wing  Loadings  (Max  of  S)  #/sq.  ft. 

NAR  ■ No.  of  Aspect  Ratios  to  be  Cycled,  Maximum  of  5 
AR(I)  - Aspect  Ratios  (Max  of  5) 

NSWP  » No.  of  Quarter  Chord  Sweep  Angles  to  be  Cycled,  Maximum  of  5 
SWP(I)  - Sweep  Angles  (Quarter  Chord)  in  Degrees  (Max  of  5) 

NTROOT  - No.  of  Root  Thickness  Ratios  to  be  Cycled,  Maximum  of  5 
TROOT(I)  - Root  Thickness  Ratios  (Max  of  5) 


TAPER  - Wing  Taper  RaCio 

SFVET  - Fuselage  Wetted  Area 

SHT  - Horizontal  Tall  Voliae  Coefficient 

LFUS  - Fuselage  Length 

CWRT  • Root  Chord  of  the  Wing  in  Inches 

MAC  - Wing  Mean  Geometric  Chord  in  Feet 

SWWET  - Wing  Wetted  Area 


$PROP 

IPROP  - 1 for  Turbofan  (P&W  401} 

• 2 for  Turbojet  (GE  J79) 

■ 3 for  Reciprocating  (Incoming  I GO-540) 

AB  - 1 for  afterburner 

AB  4 1 for  no  afterburner 

SLTH  - Thrust  to  Weight  Ratio  of  the  Vehicle 

TOWE  - Thrust  to  Weight  Ratio  of  the  Engine 

SFCF  - Specific  Fuel  Consumption  Modification  Factor 

To  end  the  input  and  begin  the  mission  simulation,  an  *E'  is 
placed  in  column  2. 

The  input  used  in  this  study  vast 
THIS  IS  THE  FIRST  TEST  CASE  FOR  ADCA 
$ DESIGN 

STOREX  - 2.,  IPRINT  - 0.,  LF  - 9.75, 

AMMAX  - 2.,  VMAX  - 794., 

INITER  - SO. , 


$END  OF  DESIGN 
$ WEIGHTS 
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WAVUN  - 1350.,  WSTORX  - 5500.,  WMA  - 600., 

GEARF  - .97,  WINGF  - .83,  BOQYF  • .85,  TAILF  - .83, 

$END  OF  WEIGHTS 
$GE0H 

NWOS  - 5,  WOSIN  - 70. ,80. ,90., 100., 110., 

WAR  - 5,  ARIN  - 2.5,  3.0,  3.5,  4.0,  4.5, 

NTROOT  - 1,  TROOTIN  - .04 

NSWP  - 1,  SWPIN  - 52.4,  TAPER  - .25,  LFUS  - 61.75, 

$END  OF  GEOM 
$PROP 

IPROP  - 1,  AB  - 1.,  SLTH  - .70,  TOWE  - 8.,  SFCF  - 1.2, 

$END  OF  PROP 
$MISSI0N 

N - 11,  GEE (7)  - 6.5,  NTURNS(7)  . .8,  OELM  - .4,  TIME(l)  - .75, 

TIME(ll)  - 20.,  Q«1N  - 2.0, 

ALT(l)  - 0.,  40000.,  40000.,  40000.,  50000.,  50000., 

50000.,  50000.,  55000.,  45000.,  0., 

AMACH(l)  - 0.,  .8,  .8,  .8,  1.6,  1.6,  1.6,  1.6,  1.6,  .8,  .4, 

RAWGE(l)  - 0.,  1.E6,  200.,  -3.,  1.1E6,  150.,  -7.,  -1.,  200.,  200.,  -5., 
$END  OF  MISSION 
E 

Sbissat. 

Th*  value*  of  th*  design  variables  and  the  corresponding  performance 
values  were  output  to  Tap*  2 and  each  value  was  of  the  format  F10.2. 

It  was  also  possible  to  have  a mission  summary  printed  out  which 
included  an  analysis  of  each  leg  in  th*  mission  along  with  a weight 


... 

lem  ii— mi>  .1  idmm 
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summary.  The  component  weights  are  also  output  but  because  of  the  many 
approximations  utilised  to  derive  these  weights,  only  the  weight  empty 
and  TOGW  values  can  be  considered  statistically  valid  (Ref  4:7).  This 
summary  data  is  available  on  Tape  6. 


Appendix  C 


degression  Analysis  Computer  Program  Input/Output 


Input 

All  numeric  data  requires  a decimal  point  with  the  exception  of 
card  number  three. 


Card  No.  Col.  No. 


Description 


1*6  The  word  PROGRM  must  appear;  this  Is  a required 


12-17  A six  character  Identification  word  In  alpha- 
numeric format 

22-31  Number  of  input  variables  per  data  point  including 
dependent  variables 

32-41  Number  of  variables  to  be  added  by  synthesis 

42-31  Number  of  variable  synthesis  cards 

52-61  Number  of  labeled  variables 

62-71  Number  of  data  points 

1-6  The  word  PRGOPT  must  appear;  this  is  a required 


12-21  Number  of  Format  cards  If  the  default  Format  of 
(11X,6F10.0)  Is  not  used 

22-25  TRUE  If  Input  data  is  on  magnetic  tape 
otherwise  leave  blank 

32-35  TRUE  rewinds  auxiliary  input  data  tape 
otherwise  leave  blank 


52-55 


62-65 


3 1-3 


a*  needed  1-6 


TRUE  prints  covariance  and  correlation  natrlces 
otherwise  leave  blank 

TRUE  calculates  the  equation  based  on  a curve 
through  the  origin  (zero  intercept) 
otherwise  leave  blank 

If  variable  synthesis  data  is  on  cards  then  the 
integer  0 must  appear;  if  the  data  is  assigned 
to  Tape  8 then  a 1 aust  appear;  this  is  a 
required  card;  the  number  must  be  in  colunn  3 
The  word  VARSYN  aust  appear;  these  are  the 
variable  synthesis  cards  which  fora  the  second 
order  teras  in  the  equations;  in  general  these 
cards  are  optional  but  were  required  in  this  study 


8-10 


Synthesis  open 

it Ion  code  which  can  take  on  the 

following  valu< 

>s: 

1.0 

VAR(I) 

Constant 

2.0 

VAR(I) 

VAR(J)  , 

3.0 

VAR(I) 

(VAR(J)) “ 

4.0 

VAR(I) 

EXP (VAR (J)) 

5.0 

VAR (I) 

SIN(VAR(J)) 

6.0 

VAR(I) 

COS(VAR(J)) 

7.0 

VAR(I) 

TAN (VAR (J)) 

8.0 

VAR(l) 

Natural  Log  (VAR(J)) 

9.0 

VAR(I) 

ARCS IN (VAR ( J ) ) 

10. 

VAR(l) 

ARCCOS (VAR(J)) 

11. 

VAR(I) 

VAR(J)  ♦ VAR(K) 

12. 

VAR(I) 

VAR(J)  - VAR(K) 

13. 

VAR(I) 

VAR(J)  * VAR(K) 

14. 

VAR(I) 

VAR (J) /VAR (K) 

15. 

VAR(I) 

VAR(J)**VAR(K)  i.e. 
VAR(J)  raised  to  the 
power  - VAR (K) 

16. 

VAR(I) 

m 

ARCTAN(VAR(J)/VAR(K)) 
(VAR(J)**2  ♦ VAR(K)**2)* 

17. 

VAR(l) 

m 

18. 

VAR(I) 

m 

MAX (VAR (J),  VAR(K),  Constant) 

19. 

VAR(l) 

m 

MIN(VAR(J),  VAR(K),  Constant) 
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12-21 

22-31 

32-41 

42-61 


as  needed  1-6 


12-15 

16-21 

22-25 

26-31 

32-35 

36-41 


Assigns  a value  to  I in  the  above  operations 

Assigns  a value  to  J in  the  above  operations 

Assigns  a value  to  K in  the  above  operations 

Assigns  a value  to  the  Constant  in  the  above 

operations 

The  word  LABELS  must  appear;  is  used  to  assign 
alphanumeric  titles  to  each  of  the  variables, 
both  input  and  synthesised;  is  an  optional  item 
in  the  program 

Index  number  for  the  first  labeled  variable 
A six  character  label  for  the  first  labeled 
variable 

Index  number  for  the  second  labeled  variable 
A six  character  label  for  the  second  labeled 
variable 

Index  number  for  the  third  labeled  variable 
A six  character  label  for  the  third  labeled 
variable 


62-65 

66-71 


as  needed  1-6 


12-71 


Index  number  for  the  sixth  labeled  variable 
A six  character  label  for  the  sixth  labeled 
variable 

The  word  FORMAT  must  appear  if  the  input  data 
points  are  in  a Format  other  than  the  default 
value  of  (11X,  6F10.0) 

Valid  Format  specification  for  the  input  data 

points;  must  begin  with  a left  parenthesis  and 

<s 

end  with  a right  parenthesis 
36 


(Data  points  ara  Input  at  this  point) 
as  neaded  1*6  The  word  PROBLM  oust  appear*  this  is  a required 

card  and  initiates  the  analysis  for  the  selected 
independent  variable 

12-21  The  index  nuaber  of  the  variable  chosen  to  be  the 
dependent  variable 

22*31  Number  of  input  or  synthesised  variables  not 
allowed  in  the  regression  equation 

32-41  Limit  the  number  of  steps  for  calculation 

42-43  TRUE  deletes  the  detailed  printout  of  steps 
otherwise  leave  blank 

32-35  TRUE  deletes  the  summary  printout 
otherwise  leave  blank 

62-65  TRUE  deletes  the  residual  printout 
otherwise  leave  blank 

as  needed  1-6  The  word  DELETE  must  appear  if  this  optional 

routine  is  used;  is  used  in  conjunction 
with  the  PROBLM  card  columns  22-31 
12-21  The  index  number  of  the  first  deleted  variable 

22-31  The  index  number  of  the  second  deleted  variable 


62-71  The  index  number  of  the  sixth  deleted  variable 

(Additional  PROBLM  cards  can  be  Input  at  this  time  to  restart  the 
regression  analysis  using  a new  dependent  variable  or  deleting 
different  variables  from  the  analysis) 

Last  1-6  The  word  FINISH  must  appear;  this  is  a required 

card  which  causes  proper  termination  of  the 

program 
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The  only  variable  synthes la  operation  code  used  In  this  study  was 
code  13. , VAR (I)  - VAR(J)*VAR(K).  This  was  used  to  form  the  squared 
and  cross  product  terns. 

The  Input  for  this  study  was  set  up  using  seven  input  variables, 
the  four  independent  design  variables  (VAR(l)  - WOS,  VAR (2)  - AR, 

VAR(3)  «•  TW,  VAR (4)  - RN)  and  three  performance  function  parameters 
(VAR (5)  - TOGW,  VAR (6)  - DISTL  - landing  distance  - , VAR(7)  - D1ST0  - 
takeoff  distance).  Regression  analysis  was  done  on  all  three  perform- 
ance functions.  See  Pig.  for  the  input.  Note  that  spacing  is  to  be 
as  described  above. 

Output 

The  output  from  the  regression  analysis  will  depend  on  the  options 
selected  in  the  input.  A complete  listing  of  the  means  and  standard 
deviations  for  the  input  parameters  is  possible  and  a printout  contain- 
ing the  covariance  and  correlation  matrices  may  also  be  obtained.  A 
detailed  printout  of  the  variables  entered  or  deleted  at  each  step  of 
the  regression  analysis  along  with  regression  statistics  will  be  re- 
ceived unless  otherwise  specified.  A point  by  point  comparison  of  the 
actual  data  with  the  calculated  data  will  be  made  and  yields  informa- 
tion indicating  the  percent  error  at  each  point  (residual  information). 
The  coefficients  for  the  selected  equation  are  printed  and  also  appear 
as  punched  cards.  The  coefficients  are  available  on  Tape  7 with  the 
first  record  having  a Format  of  (16X,  4E16.8)  and  the  following  records 
have  a Format  of  (3E16.8). 
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Appendix  D 


Optimisation  Program  Input/Output 


Input 

Card  No.  Col.  No. 
1 1-5 

6-10 

11-15 


as  needed  1-80 


as  needed  1-5 

7-8 

10-30 

as  needed  1-80 


Description 

The  number  of  independent  design  variables 
(Format  15) 

The  number  of  performance  equations  to  be  used 
(Format  15) 

The  number  of  constraint  equations  to  be  used 
(Format  15) 

The  coefficients  for  each  performance  equation 
are  input  in  the  same  format  as  they  were  output 
from  the  regression  analysis.  The  first  set  of 
coefficients  must  be  the  function  to  be  minimised. 
The  second  set  of  coefficients  must  be  the  first 
inequality  constraint  function.  When  all  the 
constraints  have  been  read  in,  the  remaining 
functions  are  input. 

Alphanumeric  character  for  the  label  of  the 
inequality  constraint  (Format  A5) 

The  inequality  characters  LE  or  GE 
The  constraint  value  (Format  G20.8) 

The  initial  values  of  the  design  variables; 
provides  the  starting  point  for  the  search 
routine  (Format  G20.8) 
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Output 


The  output  will  Indicate  the  values  of  the  design  variables  and 
the  performance  functions  at  the  minimum  point.  All  numeric  output 
is  of  the  Format  G16.8  and  is  on  Tape  6. 
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